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Abstract 

We consider the Cauchy problem posed in the whole space for the following nonlocal heat equation: 
u t — J *u — u, where J is a symmetric continuous probability density. Depending on the tail of J, we 
give a rather complete picture of the problem in optimal classes of data by: (i) estimating the initial 
trace of (possibly unbounded) solutions; (ii) showing existence and uniqueness results in a suitable 
class; (Hi) giving explicit unbounded polynomial solutions. 
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1 Introduction 

This paper is concerned with the Cauchy problem in optimal classes of data for the following nonlocal 
version of the heat equation: 

u t = J*u-u, (x, t) £ R N x R + . (1.1) 

Here, J : R N — > R is a symmetric continuous probability density and / * g stands for the convolution of 
functions / and g. Throughout the paper, the initial data 

u(x,0) = u (x), x<ER N , (1.2) 

are only assumed to be locally bounded (and not necessarily nonnegative). 

This equation and some of its variants have been studied by a number of authors recently and in 
various directions (see for instance [5] and the references therein). These works are essentially dealing 
with the class of either bounded or L 1 initial data. On the contrary, the aim of this paper is to study the 
Cauchy problem in classes of data that are not necessarily bounded, nor integrable, so that the "usual" 
tools (Fourier transform, fixed-point theorems) do not work. 

Let us mention first that in the case of continuous and compactly supported initial data uq, the 
solution of (|l.l[) - (|1.2p can be written as: 

u(x, t) = e ~ t u (x) + (u>(t) * u Q )(x) , 
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where u> is the regular part of the nonlocal heat kernel (see [5] and an expansion of uj in Section [2]). 

Now, as is well-known, in the case of the "classical" heat equation Ut = Am in R N x R+, the optimal 
class consists of initial data of functions uq satisfying the estimate: 

\u (x)\ ^ c e a ^ 2 , for some a, c > . (1.3) 

Uniqueness also holds in the class of solutions satisfying this type of estimate - see for instance [7J . In this 
paper, we shall prove that nonnegative weak solutions have an initial trace which is a locally integrablc 
function (see more precise formulations in Propositions 12.5]) . 

We would like to mention that after this paper was finished in its preprint version, we came across a 
recent article by Nathael Alibeau and Cyril Imbcrt 1 , so that both articles were written independently. 
In pQ , the authors also deal with unbounded solutions of a more general equation. They also consider other 
interesting aspects like regularization by the equation, that we do not investigate here. On the contrary, 
they essentially treat the case of the Fractional Laplacian while we are interested here in treating more 
general situations than power-type kernels. 

In order to explain more in detail the main differences between [T] and the present article, let us make 
several important remarks: 

(i) Let us first mention that in this paper, we could treat singular Lvy measures \x in the equation: 

ut{x 7 t) = / {u(x + z, t) — u(x,t)\ dfi(z) . 

Indeed, the question of obtaining the optimal behaviour of initial data is only linked to the tail of 
fj, and not to the presence of singularities at the origin. Following the notation in [T] , we can always 
decompose the Lvy measure \x as the sum of its singular part near the origin, and the bounded 
measure fib — /ul-flx^i} which conveys the tail information. In this paper, we only deal with the /im- 
part since the singular part does not play any role here. We refer to Section 17721 for some discussion 
on that issue. 

(ii) We would like to put emphasis on the fact that in the problem we are facing, one can divide the 
measures fj, (or kernels J) into two categories, according to the speed of decay at infinity. Following 
[2J|3], we say that the kernel has a slow decay if J(y) = c^^ v ^^ with u(y) -> as \y\ — » oo. 
On the contrary, the kernel has a fast decay if w(y) — > oo. The limit case of exponential decay is 
somewhat in between, but here it may be considered as being "slow decay" -type. 

(Hi) In the case of slow decay kernels, using elliptic barriers works quite well: if one can find a function 
/ ^ and a A > such that J * / — / ^ A/, then the function 

1>(x,t) := e At (l + /(a:)) 

is a natural supersolution of the problem. This is essentially what is done in [T] in the case J(y) = 
\y\~ N ~ a , for which f(x) = \x\& for any < (3 < a gives rise to a supersolution. Hence the problem 
is at least solvable in the class of initial data such that |wo(a ; )| ^ c o(l + f( x ))- We construct here 
other explicit barriers for different kernels, for instance the case of tempered a-stable laws (and 
CGMY-processes, see [4]) for which J(y) — e ~\ z \ /\z\ N+a , see Theorem II .31 

(iv) A natural related question is to prove that the barriers actually give the optimal behaviour for the 
initial data. Although this is not stated in [T] as such, the existence result obtained there is optimal 
for power- type kernels (see Section 14.11 below). In all the slow decaying kernels we consider, we 
prove that the explicit elliptic barriers we construct are optimal. 
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(v) On the contrary, in the case of fast decaying kernels the elliptic barriers do not give the optimal 
class of initial data, which does not allow to treat those cases in full generality. Instead, one has 
to use another approach by estimating directly the non local heat kernel in a different way and 
construct other supersolutions than the ones in [Tj, more adapted to the problem. We consider, the 
typical cases of gaussian laws and compactly supported kernels. 

(vi) We also investigate related problems like the phenomenon of blow-up in finite time, asymptotic 
behaviour of unbounded solutions and we construct explicit polynomial solutions. 

Here are now our main results. 

Fast decaying kernels - Let us take a symmetric kernel J such that the support of J is exactly 
Bi = {\x\ ^ 1}. It is known that in this case, equation ([l.ip shares a lot of properties with the local heat 
equation, [5]. So, is the class of optimal data for also given by fj 1 . 3[) ? The answer is no: instead 

of the quadratic exponential growth we get a |a;| In |x|-exponential growth, but the result is similar in its 
form: 

Theorem 1.1 Let J be a radially symmetric kernel, supported in B\, let uq be a locally bounded initial 
data and Cq > arbitrary. Then the following holds: 

(i) If \uo(x)\ ^ coe ' 1 ' 1 "' 1 ' for some < a < 1, then there exists a global solution of (jl . 1|) — (ll.2[) . 

(ii) Ifuo(x) > coe^l 1 "^ for some fJ > 1, then there exists no nonnegative solution of (|1.2p . 

(Hi) For any < a < 1, there exists at most one solution such that \u(x,t)\ ^ C^e ^' 1 "' 1 ', where C(-) 
is locally bounded in [0, oo). 

Let us mention that this theorem cannot be obtained by using elliptic barriers, since they do not give 
the optimal (|a;| In |a;|)-behaviour. At least, for any a > 0, the function e"' 1 ' is an elliptic barrier which 
is enough to get uniqueness results in the class of solutions and initial data that grow at most like 
exponentials. On the other hand, trying a barrier like f(x) — e a \ x \ ln \ x \ leads at best to the estimate 
J * / — / ^ A|x| Q /, which does not allow to construct a suitable supersolution ip as above. 

Thus, this theorem comes from a direct estimate of the "heat kernel" u>(x,t), using its semi-explicit 
series expansion (|2.2[) . see Section [5] 

Blow-up in finite time - In the case of the usual heat equation, it is well-known that if the initial 
data has a critical growth of the order of a quadratic exponential, then the corresponding solution blows 
up in finite time. We prove here that it is also the case for gaussian kernels: 

Theorem 1.2 Let J be a centered gaussian with variance a 2 ~ 1/2, let uq be a locally bounded initial 
data and Co > arbitrary such that 

u (x) = c e WP" W) 1/2 +/(W) for some a < < /?, as ^ f(s) ^ (3s . 
then the minimal solution of (|l.l[) - ([1.2p blows up in finite time. 

For more complete results on existence, nonexistence and uniqueness, see Theorem 15.111 We conjecture 
that Theorem 1 1.2 1 is also valid for compactly supported kernels with a (|x| In |a:|)-critical rate instead (and 
more general kernels, if not all of them), we hope to address this question in a future work. 

Optimal trace for slow decaying kernels - We prove that the estimates given by the elliptic 
barriers are optimal in a number of slow decaying cases. Let us just illustrate this with a result concerning 
tempered a-stable law-type kernels: 



3 



Theorem 1.3 Let J be a kernel satisfying 

a = sup | a > : / J(y) e M (1 + \y\) N+a dy < oo} < oo . 
Let Uq be a locally bounded initial data and Cq > 0. Then the following holds: 

(i) If \uo(x)\ ^ Co e l x '(l + |x|) Ar+Q for < a < ao, then there exists a global solution of (11.21) . 

(ii) If uq(x) ^ Coe' x l(l + |x|) Ar+a for a ^ ceo, i/ien i/iere exists no nonnegative solution of p. 21) . 
(m) For any < a < ao there exists at most one solution such that 

\u(x,t)\ <C Coe' x '(l + |x|) Ar+Q locally uniformly in [0,oo). 

This theorem applies for J(y) = e~' y '/(l + |y|) Ar+Q , but also with some adaptations to the singular case 
J(y) = e /\y\ N+a , as we mentioned at the beginning of the introduction (see also Section ITT2")) . 

Asymptotic behaviour - It is well-known that if the initial data is integrable and bounded, then u(t) 
converges to zero as t —> oo (see 5 ). But if we consider solutions given by an initial data such that 
uq(x) — > +oo as \x\ — > oo, what can we get as asymptotic behaviour? We prove here that there exist 
explicit polynomial solutions which are similar to the ones for the local heat equation. Moreover, these 
solutions tend to +oo as t — > oo. For instance: 

Theorem 1.4 Let J be a kernel with finite second-order momentum ni2 > 0. Then if uq = \x\ , the 
solution of (11.11) - (|1.2j) has the explicit form 

u(x, t) — \x\ 2 + m2t. 

More generally, if J has a momentum of order 2p and Uq = \x\ 2p , then 

u{x,t) = \x\ 2p + m 2p t p + o(t p ) as £ ->• oo . 

To end this introduction, let us mention that more general nonexistence result for changing sign 
solutions can be proved provided the negative (or positive) parts of the initial data and the solution 
are controlled, see Remark |4.4I Also, similar results are easily obtained for only locally integrable data, 
provided they satisfy the various estimates for \x\ large. Indeed, this is just a matter of decomposing uo 
as the sum of a locally bounded function and an integrable function for which we can solve the problem 
separately, using the linearity of the equation. 

Finally, let us mention that the bounds we obtain are related to the large deviations estimates given 
in [21 [3] , because in both cases we are somehow estimating the nonlocal heat kernel of the equation, see 
also Section [7] 

Organization — Section [5] shows that the initial trace exists for nonnegative solutions and give some 
basic regularity results. Then we show in Section [3] general existence and uniqueness results provided 
there exists a positive supersolution of the problem. Section|4]is devoted to apply this to various examples 
of slow decaying kernels, and prove that we obtain optimal existence and uniqueness results. In Section[5] 
we deal with fast decaying kernels. In Section [5] we construct explicit polynomial solutions which tend 
to +oo as t — > oo if the initial data is unbounded. Finally we end with a section of comments, further 
results and possible extensions. 

2 Preliminaries 

First of all let us fix some notation that will be used in what follows: 
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(i) We denote by B r be the ball of radius r centered at and by \ r its characteristic function. 

(ii) Throughout the paper, J : R N — > R is a nonnegative symmetric continuous function such that 
Jrn J(y) dy = 1; the notation J* n denotes the convolution of J with itself (n— l)-times, n counting 
the number of J's, so that by convention J*° = Sq , J* 1 = J. 

(Hi) If r : M. N — > M is a smooth, nonnegative and compactly supported in B\ function, such that J t = 1, 
then r n stands for the resolution of the identity given by r n (x) = n N T(nx). 

(iv) If /, g : M. N ->lwe write / -C g if f(x) = g(x)e(x) with e(x) — > as |x| — > oo. But moreover we 
extend it as follows for time-depending functions: 

Definition 2.1 If u,v : R N x [0, oo) — > M., we say that u -C v locally uniformly in [0, oo) if 

\u(x,t)\ ^C(t)e(x)\v(x,t)\, 
for some C € L% c ([0, oo)) and e(x) — >• as \x\ —> oo . 

This definition will be useful in the construction of a supersolution ip which will tend to +oo as t — > oo. 
Notice that definition does not require to have a uniform control for all t 0. 

2.1 The nonlocal heat kernel 

It is shown in [S] that if both the initial data and its Fourier transform are integrable, then the solution 
to (jl.ip ~ (jl.2p is unique and given by: 

u(x, t) = e ' l u (x) + (uj(t) * u ) (a;) , (2.1) 

where ui is smooth. Of course, this is true for instance if the initial data is in L 2 (R JV ) or bounded and 
compactly supported (we shall use this result in various constructions in the sequel). The following lemma 
collects some basic properties of the regular part of the nonlocal heat kernel: 

Lemma 2.2 The function lo is given by 

Moreover, lo is a solution of the following problem: 

ojt — J * uj — oj + e ~ l J , 
oj(x 7 0) = 0. 

Proof. If uq and its Fourier transform are integrable, we can write the solution in frequency variables 
as follows: 

oo n 

u(z,t)=e- t (i + Y,^J n (0)MO- 

n=l 

Hence, going back to the original variables we get 

/ °° t n \ 

u(x,t) = e-*f<5o(a;) + ^ — J* n (x)j *u (x) = e -t n (a;) + (oj(t) * u )(x) , 

n=l U ' 
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(2.2) 



(2.3) 



where uj is given by (I2.2p . Notice that since ||J*™||i = \\J\\i — 1, then the series converges in L 1 . 
Moreover, a direct computation shows that to(x,0) — (in the continuous sense) and 

/ °° nt n ~' 1 J* n \ / °° t n J*("+ 1 ) \ 

^ = -^+ c ' t ( J +E n[ )=-^+ c ' t ( J + E n\ ) 

n—2 n—1 

= —uj + c ~ f ( J + J * j — ) = ^*w-cj+e~*J. 

n=l U ' 

This ends the lemma. □ 
2.2 Weak and Strong Solutions 

Let us now specify what are the notions of solution that we use throughout the paper: 

Definition 2.3 A weak solution of (jl.ll) is a function u € Lj L oc (R Ar x R + ) such that the equation holds 
in the sense of distributions. 

We may also, as usual, consider strong and classical solutions: 
Definition 2.4 Let u G Ll oc (R N ). 

(i) A strong solution of (jl.ip is a function u £ C°([0, 00); L 1 1 oc (R JV )) such that u t , J * u £ Lj oc (M JV xl + ) , 
the equation is satisfied in the Lj oc -sense and such that u(x,0) = uq(x) almost everywhere in M . 

(ii) A classical solution of ([l.ljl is a solution such that moreover u,ut, J * u E C°(R. N x [0, 00)) and the 
equation holds in the classical sense everywhere in R N x [0,oo). 

(Hi) A sub or supersolution is defined as usual with inequalities instead of equalities in the equation and 
for the initial data. 

The following result proves that in fact, nonnegative weak solutions are strong. This will allow us 
to consider only strong solutions in the rest of the paper, even for changing sign solutions which will 
be constructed as the difference between two nonnegative solutions. But let us also mention here that 
actually, the same result will be valid for changing sign solutions provided they belong to the uniqueness 
class, see Section [7] 

Proposition 2.5 Let u be a nonnegative weak solution of (jl.ll) . Then u has an initial trace u(x,0 + ) 
which is a nonnegative Lj^ function, and u is a strong solution of (11.11) . If moreover u(x,0 + ) is contin- 
uous, then u is in fact a classical solution. 

Proof. We consider the auxiliary function v(x, t) := e t u(x,t) which satisfies: 

d t v(x, t) = J * v ^ . 

This proves that the limit v(x,0 + ) is defined. Now, since u G Lj oc (R N x R+), then u(-,t) € L 1 1 oc (R JV ) for 
almost any t > 0, so that for such a time, 

< u(x, 0+) = v(x, 0+) < v(x, t) = e l u{x, t) e l^ oc (R N ) . 

This proves that the initial trace of u is indeed a Lj oc function. Even more, this also proves that 
u(-,t) 6 L 1 1 oc (R Ar ) functions for all t > 0. Now, we can write 

v(x,t) =v(x,0 + ) + [ (J * v)(x, s) ds , (2.4) 
Jo 
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which implies that J * v and then J * u are L 1 1 oc (R Ar x R + ) functions. Since u is already a Lj" oc function, 
then ut is also a Ll oc function. Hence it is a strong solution. Finally, if u(x, + ) = v(x, + ) is continuous, 
then clearly by (I2.4[) . v is continuous in space and time, up to t — 0. Hence u is continuous and the 
equation holds in the classical sense. □ 

The following technical trick will be useful in order to get some continuity for comparison results: 

Lemma 2.6 Let u be a strong solution to (|l.ip - ([1.2[) . Then for any smooth and compactly supported 
function r : R — > R, the function 

u T (x,t) :— (t * u)(x,t) — / t(x — y)u(y,t) dy 

Jr n 

is a classical solution of (|1.1|) - (|1.2I) with initial data r * Uq. 

Proof. Since r is compactly supported and smooth, the convolutions r * u$ and u T are well-defined and 
moreover, elementary properties of the convolution show that 

dt(u T ) = t * u t = t * (J * u) — t*u = J*u t — u T . 

Finally, u T is clearly continuous hence it is a classical solution and its initial trace is r * uq . □ 

Of course, the same result is valid for sub/super solutions: since r„ is nonnegative, the convolution 
maintains the inequality. 

2.3 Supersolutions 

In order to establish a full theory of existence and uniqueness for solution of (|l.ip - (|1.2p we will use special 
supersolutions tp that verify: 

i/> G C°(R W x [0,oo)),V > 0, 

ip(x, t) — > +oo as |i| — > oo uniformly for t £ [0, oo) , (2-5) 
rp t > J * ip - ip ■ 

A typical way to construct such supersolutions is to use barriers, i.e., functions / : — > R continuous, 
nonnegative, with f(x) — > oo as \x\ — > oo and such that for some A > 0, 

■/*/-/< A/, A>0. (2.6) 

Then ip(x,t) := e xt f(x) satisfies (|2.5I) . For instance, if J has a finite second-order momentum m2, then 
/(x) = x 2 fulfills the requirements, with A = 7712, this is a calculus similar to the one performed in 
Section [5] 

Another approach that we shall use in the case of fast decaying kernels is the following: since ui 
satisfies (|2.3I) . then it is a supersolution of equation (|l.ip . Thus, if / is nonnegative and such that the 
spatial convolution of to and / converges, then the function (x, t) h-> (ui(t+l)*f)(x) is also a supersolution. 

3 Existence and uniqueness of unbounded solutions 

3.1 Comparison 

Comparison in the class of bounded solutions is well-known if J is compactly supported. We extend now 
the comparison result to more general kernels and not only for the class of bounded solutions. To do so, 
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we need to control the points where the maximum is attained, which is done by using the supersolution 
ip defined above. 

Proposition 3.1 Letip satisfy (|2.5p . Let u be a strong subsolution of p. ID and u a strong supersolution 
of such that u(x,0) ^ ii(x, 0). Ifu—u <C ip locally uniformly in [0,oo) thenu ^ u almost everywhere 

in R N x R + . 

Proof. We shall do the proof in two steps: we assume first that the sub and supersolutions are continuous, 
and then use a regularization procedure to extend the result. 

Step 1 - Assume that both u and u are continuous and consider the function 

w s (x,t) — e t (u(x,t) — u(x,t) — 8ip(x,t) — 8) , 
with 8 > 0. This function satisfies 

d t {w s ) -J*w s ^0. 

Let (xo,to) be the first point at which w s reaches the level —8/2 and assume that to is finite. Since 
u — u<S.ip locally uniformly, we get that xo < oo. Moreover, since the function w s is continuous and 
w 5 (x,0) sC -8, then t > 0. 

Notice that J * w 5 ^ —8/2 for t 6 (0, to). Then, dt{w 5 ) ^ for t € (0, to), which is a contradiction. 
Therefore, w s $S —5 for all time t > 0. Finally taking <5 — > 0, we obtain the desired result. 

Step 2 - Now we need to relax the continuity assumption. To this end, we use Lemma |2~51 with a suitable 
resolution of the identity r„ (see preliminaries) to see that u n :— r„ * u is a continuous subsolution and 
similarly, u n :— r n *u is a continuous supersolution. Moreover, since r n ^ 0, their initial data are ordered: 

U n (x, + ) = r n * u < t„ * m = w„(x, + ) . 

In order to use Step 1 for the convolutions u n and U n , let us notice first that if i\> satisfies (I2.5p . then 
:= T n * ip is also a supersolution (again, this is just convoluting the inequation). We only have to 
verify that ip n locally uniformly. 

Actually, this comes from the fact that by definition, u — u <^ ip and since we take the convolution 
with a nonnegative r„, 

u n - u n < c(t)r„ * {e(-)ip(-,t))(x) < e'(a;)(T n * ip)(x) , 

where e'(x) = max{e(y) : \y — x\ ^ 1/"} still verifies e'(x) — > as |x| — > oo. Hence, u n — u n -C an d 
we apply Step 1 with u n , u n and i/in, which yields: 

u n ^ ?2 n in R N x [0, oo) . 

Finally we pass to the limit as n —¥ oo. Since u and u are locally integrable, we get u ^ u almost 
everywhere. □ 

As direct corollary, we have the following uniqueness result, which is valid only in a suitable class of 
solutions, as it is the case for the local heat equation. 

Theorem 3.2 Let ip satisfy (|2.5[) . Then there exists at most one strong solution u of (|1.2[) such 

that \u\ -C ip locally uniformly in [0,oo). 

Proof. If two solutions u and v satisfy |u|, |u| <C ip locally uniformly in [0,oo) then we have at the 
same time u — v <Si ip and v — u <C ip locally uniformly in [0, oo). Using these solutions as sub- and 
supersolutions allows us to obtain u ^ v and v ^ u so that they are identical. □ 
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3.2 Construction of a solution 



Theorem 3.3 Let %p satisfy (|2.5p . Let uq be locally bounded and assume that for some cq > 0, 

\u (x)\ < coip(x,0) ml". 
Then there exists a solution u of (ll.ll) - (|1.2l) . which is given by 

u(x, t) = e ~*uo(x) + (w(t) * Uo)(x). 
Moreover, the following estimate holds: \u(x,t)\ ^ coip(x,t) in M. N x R + . 

Proof. We separate the positive and negative parts of the initial data (ito)+ and («o)-, an d solve 
separately the two problems by linearity of equation Hence, let it+ be the unique solution, see [5], 

of the truncated problem 

(u n ) t = J *u n — u n , (x,t) G R™ x (0,oo), 

u n (x,0) = (u ) + (x)xn(x), x e R". 

Since, u n (x, 0) € L 2 (R Ar ) the representation formula (|2.ip holds. Moreover, since M n (ir, 0) € L°°(R Ar ), 
then 

U+(M) = e-*u n (x,0) + (w(t)*w n (:r,0))(x) s: ||u n (a;, O)^ . (3.1) 

So, we have obviously u+ <C ip locally uniformly in [0,cx)). Therefore, we can apply the comparison 
principle, Proposition 13. 1[ to get 

(i) the sequence {u^} is monotone nondecreasing; 

(ii) u+(x,t) < c Q ip(x,t). 

Hence, there exists a limit u + defined in all x (0,oo). We have to check that it is in fact a 
solution of (11.11) . Passage to the limit in the equation is done by using the dominated convergence for 
the convolution term: the limit function u + is in L 1 1 oc (R Ar x (0, oo)) and J * it+ converges to J * u + 
in L 1 1 oc (R A '), thus we recover a strong solution with initial data (ito)+ . Moreover, passing to the limit 
in (|3.1[) and in point (ii) above we obtain 

u + (x,t) = e~ t (u ) + (x) + (w(t) * (u Q ) + )(x) , u + (x,t) ^ c ip(x,t) . 

The same construction, using again cqi/j as a nonnegative supersolution gives a solution u~ with initial 
data (ito)-, the negative part of u$. Thus, u = u + — u~ is a solution with initial data uq = (?io)+ — (ito)-- 
Comparison with ip is straightforward. □ 

Remark 3.4 Several remarks are to be made now: 

(i) The solution u + constructed in this theorem can be used as a minimal solution in the class of 
nonnegative solutions. Indeed, assume that v is another nonnegative solution with the same initial 
data Mo ^ 0. Hence v is also a supersolution and it can be used to play the role of cgip- Thus, 
u+ ^ v and passing to the limit we get u ^ v. 

(ii) There is a gap between the uniqueness and existence theorems: one requires that u <C ip to get 
uniqueness while we are only able to construct a solution comparable to ip. Thus, if one wants a 
global result of existence and uniqueness, one shall use two functions ipi <§; ip2 satisfying both (|2.5[) : 
there exists u ^ ipi and it is unique, because u ^ ipx <C ip2- 



9 



We end this section with a first estimate of the initial trace, using the same construction as for the 
existence theorem. Notice that this estimate only concerns nonnegative (or nonpositive) initial data, but 
since \uq\ = — Uq , then we shall be able to use the argument to obtain some precise estimates also 
for changing sign solutions in the next section. 

Corollary 3.5 Let u be a nonnegative solution of (|1.1[) - (|1.2[) . Then the initial trace uq = w(-,0 + ) ^ 
of u satisfies the estimate: 

For any n ^ 1, (J* n * uq)(x) < oo . 

Proof. We use the same construction as above, considering UoXn as initial data. The sequence u n of 
solutions with initial data UoXn is monotone nondecreasing (recall that here uq ^ 0). Moreover, since 
UoXn € L 2 (M Ar ) the representation formula holds: 

u n (x,t) = e~ t u (x)xn{x) + (uj(t) * (u Xn))(x). 

Since there exists a solution u, then the minimal solution u constructed as in the previous proof, as the 
limit of u n , has to be finite almost everywhere. This implies that 

(uj(t) * wo) (x) < oo a.e. in R w , 

and actually everywhere since uj is continuous. Using the explicit formula for w, see (|2.2I) . we deduce that 
for any n ^ 1, we have necessarily J* n * uq < oo everywhere in M. N . □ 

4 Optimal behaviour of the initial data for slow decaying kernels 

Using the results of the previous section we are able to find the optimal class of initial data such that 
there exists a unique solution to (| 1 . 1 [) — (jl.2l) . As we pointed out in the Introduction, the main point 
here is that J decays slow at infinity. In fact, we will consider three different J's, all of them decaying 
at infinity slower than exp(— a > 0. All three cases are treated in the same way: using elliptic 
barriers. That is, the main point consists in finding the biggest function / that verifies hypothesis (|2.6[) . 
The existence of such a function implies the existence of the supersolution ip. 

4.1 Power-type kernels 

As mentioned in the Introduction, we recover in this subsection the existence and uniqueness results 
in [I] for power-type kernels. We also give complementary results concerning nonexistence of solutions 
and optimality of the barriers. 

Let us suppose that 

70 := sup 1 7 > : J J(y)\yp dy < oo j e (0, oo) . (4.1) 

That is, let us assume that J has momentum of order 7 for all 7 < 70- Observe that, by continuity of 
the function 7 h> f J{y)\y\ 1 dy, we get that the momentum of order 70 is infinite. Kernels satisfying 
condition (|4.1j) are for example those with potential decay J(x) ~ | cc | ("yo-t-^v) f or j^j l ar g e . 

Lemma 4.1 Let J be a kernel satisfying (|4. 1 [) . For any 7 < 70 there exists A = A(J, 7) such that J 
verifies (|2.6|) with f(x) = 1 + |x| 7 . 
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Proof. Since \x — y\ ^ 2max{|a;|, \y\}, we have 



J*(l + |:r| 7 )-(l + M 7 ) = J*\x\ 



1 



|.t| 



<J*|af= f J{y)\x-y\^dy 





J(y)\x - yr dy + / 



K N \B(0,|a;i) 



J{y)\x - y^dy 



^M 7 2 7 / J(y)dy + 2 



J(y)\yVdy 



JB(0,\x\) 



R N \B(0,\x\) 



< A(|zp + 1), 



which holds if A ^ 2 7 max 



J(2/)M 7 d2/}. 



□ 



Theorem 4.2 Let J be a kernel satisfying (|4.1[) . Let «o &e a locally bounded initial data such that 
\ u o( x )\ ^ c o(l + M 7 ) / or Co > and < 7 < 70. Then, there exists a global solution of (|1-1[) - (|1.2I) . 
which satisfies 

\u(x,t)\ <C Co(l + M 7 ) locally uniformly in [0, 00). 
Moreover, u is the unique solution in this class. 

Proof. It follows directly from theorems 13.21 15751 and Lemma I4~T1 by choosing V> = coe A *(l + |a;| 7 ). □ 

Theorem 4.3 Let J be a kernel satisfying (|4.ip . Assume thatuo(x) ^ co(l + |cc| 7 ) for cq > and^ ^ 70. 
Then there exists no nonnegative solution of (|l.ip (jl.2l) . 

Proof. By Corollary [531 if there exists a nonnegative solution, then J *uq < 00, which is a contradiction 
with the growth of uq. □ 

Remark 4.4 If uo is a changing sign initial data, the non-existence theorem holds in the class of solutions, 
whose negative (or positive) part is somehow controlled. More precisely, if we consider initial data uq 
such that J * (uo)_ < 00 and J * (uq) + — 00, then there are no solutions in the class {w_ ^ ip}. Indeed, 
assuming there exists a solution u in this class, then we can consider the nonnegative supersolution 
v = u + tp, which allows us to construct the minimal solution u with initial data vq = uq + ^(0) as in 
Theorem 13.31 But then, J * vq has to be finite, which contradicts the hypothesis on the initial data since 
of course J * ip(0) < 00. 

4.2 Exponential-type kernels 

If the initial data is an exponential type one, the analogous of having momentum of order 7, see (|4.1I) . is 



For instance, this is the case for J{x) ~ coe T°l x l. 

Lemma 4.5 Let J be a kernel satisfying (|4.2I) . For any 7 < 70 there exists A = A(J, 7) such that J 
verifies (|2.6[) with f(x) = e 7 l x l. 




(4.2) 
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Proof. The proof follows as before. We only note that \x — y\ — \x\ ^ \y\, then 

J* e 71 * 1 - e 7 ^ 1 ^ J* e 7 ^ = e 7 ^ / J(y) e^ x - y ^^ dy 

< e 7|a;| / J(y)e^dy 
^ Ae 7 '* 1 , 

which is true, if A > J RN J(y)e J ^ dy. 

Hence, following the same arguments given in the previous subsection, with ip(x,t) 
obtain the following results: 

Theorem 4.6 Let J be a kernel satisfying (|4.2[) . Let uo be a locally bounded initial data such that 
\ u o(x)\ ^ coe 7 ' 2 ' for cq > and < 7 < 70. Then, there exists a global solution of p.lj) - (|1.2p . which 
satisfies 

\u(x,t) \ <C coe 7 ' 2 ' locally uniformly in [0,oo). 
Moreover, u is the unique solution in this class. 

Theorem 4.7 Let J be a kernel satisfying (|4.2[) . Assume that uq(x) ^ coe 7 '^' for cq > and 7 ^ 70- 
Then there exists no nonnegative solution of (|1.2I) . 

Moreover, from Remark 14.41 we also have nonexistence for solutions with changing sign if we control 
the negative part of the initial data. 

4.3 Tempered a-stable type kernels 

Assume now that the kernel J verifies 

70 = sup { 7 > : J J(y) e 7|y| dy < c^} e (0,oo), (4.3) 
a = supja > : / J(y)e 7o|y| (l + \y\) N+a dy < 00} G (0,oo). (4.4) 

e -to\x\ 

This is the case for instance if we consider J (a;) ~ ; — rrrr; — as \x\ — > 00. 

w (1 + \x\) N+a o 1 1 

It is easy to check that f(x) — e 7 ^', with 7 < 70 is a barrier, but it is not optimal as the following 
lemma shows: 

Lemma 4.8 Let J be a kernel satisfying (|4.3|) . For any a < ao there exists A = A(J, a) such that J 
verifies (JUD with f(x) = e 7n l x l(l + \x\) N+a . 

The proof follows the same lines as in Subsection 14.11 We prefer to write here the power as N + a 
instead of a, in order to be more consistent with the usual notation for tempered a-stable laws. 

Proof. Let us check that f(x) = e l0 ^(l + |x|) A ' +Q is in fact a barrier: 

J* e"*>W(l + \z\) N+a - e 7oN (l + \x\) N+a ^Xe^ lxl (l + \x\) N+a , 



□ 

= e^e^ 1 ! we 
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for some A to be chosen later. Indeed, Since / is a positive function and \x — y\ — \x\ $J \y\ we have that 
J*f-f<J*f= [ J(y)e^-y\(l + \x-y\) N+a dy 

= e 7o|a;| / J(y)e^ x ' v ^^{l + \x-y\) N+a dy 

JR™ 

< e^^l / J(y)e^(l + \x~y\) N+a dy=e^I. 

JR™ 

In order to treat the power term, we note that 1 ^ (1 + |a; — y\) ^ 2(1 + niax{|cc|, \y\}), then 
/= / J(y)e^(l + \x-y\) N+a dy + [ J(y) e (1 + \x - y\) N+a dy 

JB(0,\x\) JR n \B(0,\x\) 

^2 N+Q (l + \x\) N+a f J{y)e^dy + 2 N+a f J(y) e 7ol ' y| (1 + \y\) N+a dy 

Jb(0,\x\) JS, n \B(0,\x\) 

^2 N+a (l + \x\) N+a ( [ J(y)e^dy + f J(y) 1(1 + \y\) N + a dy) . 

K JB(0,\x\) JR N \B(0,\x\) 7 

Hence, using fO} it is enough to take A ^ 2^+" J RN J(y)e^(l + \y\) N+a dy to get the result. □ 
The results of existence, uniqueness and non- uniqueness are straightforward. 

Theorem 4.9 Let J be a kernel satisfying (|4.3[) . Let uq be a locally bounded initial data such that 
\ u o( x )\ ^ coe 7 '^^! + I^D^^" for cq > 0, < 7 < 70 and < a < ciq . Then, there exists a global 
solution of p.ip ~ (|1.2j) . which satisfies 

\u{x,t)\ < c e 7|:r| (l + \x\) N+a locally uniformly in [0,oo). 

Moreover, u is the unique solution in this class. 

Theorem 4.10 Let J be a kernel satisfying (|4.3[) . Assume that u (x) ^ c e^ x \(l + \x\) N+a for c > 0, 
7 ^ 7o o,nd a ^ a>o- Then there exists no nonnegative solution of (jl.ip - (|1.2p . 

5 Optimal behaviour of the initial for fast decaying kernels 

We face now the cases of fast decaying kernels. Although in these cases the constructions of barriers is 
also possible, we do not get the optimal behaviour and we need to develop a different approach. 

The main point here consists in estimating the nonlocal heat kernel, 

^ — ' n\ 

71=1 

associated to the equation, in order to get the optimal behaviour for the initial data. 
5.1 Compactly supported kernels 

In this section we assume that J is radially symmetric and compactly supported in B p for some p > 0. 
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Proposition 5.1 Let J be a kernel such that supp(J) = B p for some p > 0. Then there exists two 
constants 03,04 > and for any < a < p there exist two constant C\,C2 depending on a such that 

Cl e "* e -(V^M 1" |*| e (c 2 +lnt)|x| ^ ^ t) ^ ^ e -t e -(l/p)|*| In |x| e (c 4 +lnt)|x| _ (g^ 

We begin with the upper estimate: 

Lemma 5.2 Let J be a kernel such that supp(J) = B p for some p > 0. Then there exist c 3 , c 4 > such 
that the following upper estimate holds: 



00 j-n 
Z — ✓ *il 



n=l 

Proof. First it is clear that since J has unit mass, then ||J* n ||oo ^ Halloo = c > 0. Second, since the 
support of J* n is exactly i?( n+1 ) p , then it is enough to estimate, for fixed t > 0, the function 

j-n 

where [r] stands for the entire part of r > 0, that is, the greatest integer below (or equal to) r. Using the 
Taylor expansion of the exponential, we get that 



y-= f 



cSds ^^7T- (5-2) 



Then we use Stirling's formula for the factorial: K\ = y/2irK(K/e) K (1 + 0(1/K)) to obtain that there 
exists d > such that: 



n! 



Coming back with K = \\x\/ p\ wc then obtain: 



V" *2 < cc^IM/p] e -(V2) ln[|x|/p] e -[|x|/p](ln[|x|/p]-l) 



n^[|x|/p] 

Then it follows that for some C3, c 4 > 0, 

_ J* n ( x ) ^ C 3 e "(V/ONI ln k| + (c4+lnt)|a;|/p ^ 

n! 

n=l 

which is the desired estimate. □ 

Lemma 5.3 For any < a < p, there exist c > and fi e (0, 1) depending only on J and a such that 

Vn^l,Vx€B na , J m {x)^cp n . (5.3) 

Proof. For some technical reason, we need to introduce a' = (a + p)/2 e (a, p). Since J is symmetric 
and compactly supported in B p , noting y = (yi, . . . , y^) we set 

p:= J(y)dy>0, 

J {<?' ' <Vn <p] 
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which is invariant under rotations of the coordinates. Recall also that supp( J* n ) = B np which strictly 
contains B na , so that for no ^ 1 to be fixed later, we can define 

c := min J* n " / p > 0. 

With this choice, property (|5.3I) is valid for n = no and we shall show the property by induction for 
n > no- So we us assume that in B nc7 , J* n c/i™ and for x € Bi n+ x\ a let us estimate: 

j*(n+i) (a;) = f J( x - y )r n (y)dy= [ J(y)J* n (x-y)dy. 

We assume without loss of generality that x = (0, . . . , 0, Xn) 

J< n+1 \x) > / J(y)J*"(x-y)dy. 

{c r '<ajv<p} 

Then since |x| < (n + l)a, if y satisfies a' < yjy < p, we have: 

(N-l) 

\x- y\ 2 =\x N - vn\ 2 + iy*i 2 ' 

»=1 

<((n + l)cx-c/) + p 2 , 
^(na + a - a 1 ) 2 + p 2 . 

Since a < a' we have for n big enough that (na + a — a') 2 + p 2 ^ (na) 2 and it follows that for such x, y, 
\x — y\ ^ na. Hence we fix no = no(a, a' , p) such that for any n ^ no and x € Br n+ \\ a we have: 

J< n+1 \x) > [ (cp n )J(y) dy = cp n+1 . 

J {<?' <VN<P} 

Thus the property is proved by induction for n > no- Finally, the cases when n ^ no only concern a finite 
number of terms so that, up to redefining c and p taking the min over all the first terms, we obtain (|5 .3[) 
for any n 1. □ 

Proof of Proposition [5TT1 The upper estimate follows directly from Lemma \5. 21 so we turn to the 
lower estimate. We first notice that for any fixed n 

. . t t"J* n (x) 
n\ 

If we take n = [\x\/a] — 1 then x G B na and we may use Lemma 1531 there exists c, p > depending on 
a such that 

w(s '* )>ce "-iH7^r- 

Thus, using again Stirling's formula for the factorial, we obtain that for some c, d > 0: 

uj[x,t) > ce-^-^/^kl^lxl+fc'+H^lxl _ 

Hence, we have obtained indeed that there exists ci,C2 > such that 

w(af, t) ^ ci e -* e -(V<0M In |*| e (c+lnt)M ^ 

which ends the proof. □ 
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Theorem 5.4 Let J be a radially symmetric kernel, compactly supported in B p for some p > 0. Let Uq 
be a locally bounded initial data and Co > 0. Then the following holds: 

00 If\u (x)\ < coe^l^N f or 

some < a < 1/p, then there exists a global solution of (11.21) . 

(ii) Ifu {x) ^ coe^l^N / or 

some j3 > 1/p, i/ien there exists no nonnegative solution of (|l.ip - (|1.2j) . 

Proof. We consider the sequence {u n }, given by the solution of (jl.ip with initial data u n (x, 0) = uq(x)xu- 
This functions are given by convolution: 

u n (x, t) = e ~*u o (a:)Xn(20 + * MoXn)(a). 

(i) - If 1^0(^)1 ^ coe"' 1 ' 1 "' 1 ' for some < a < 1/p, then oj(t) * uq is defined for any t > 0, and 
the sequence is uniformly bounded by 4>(x,i) — e~ t uo(x) + uj(t) * uq. Therefore, following the 
same argument given in the proof of Theorem 13.31 we get a global minimal solution, which is given by 
convolution. 

(ii) - Let us assume that uq(x) ^ coe^' 1 ' 1 "' 1 ' for some (3 > 1/p, and that there exists a nonnegative 
solution u. Then notice first that there exists also a a € (0, p) such that (3 > 1/a and either the negative 
or the positive part of uq satisfies a similar estimate, for instance: 

(uo)+(aO ^coc^l 1 "!*!. 

We use an approximation by compactly supported initial data (uo)+x« an d we see that by comparison, 
in R N x R + we have 

e "*(w )+Xn + * (u ) + Xn < M < 00 . 

But on the other hand, since j3 > l/o, u(t) * (iio)+Xn ~* +°° as n °°! which is a contradicton with 
the bound above. The conclusion is that no solution can exist. 

□ 

Remark 5.5 Several remarks are to be made: 

(i) Actually the estimates on ui allow us to treat more delicate situations like for instance (taking p = 1 
for simplicity), uq(x) ~ e ln bl-bl ln ( ln kl) which gives rise to a global solution. 

(ii) Finer estimates (in x) can be obtained by using a better expansion of the factorial - which is known, 
but anyway, the (c + lni)|a:|-terms cannot be avoided which implies an estimate of the order of 
in the exponential. 

(Hi) An interesting question concerns the critical growth, when uq(x) behaves like e ln \ x \ . We 
shall see in the gaussian case below - see Theorem 15.111 (Hi) - that solutions with critical initial 
data blow up in finite time, as it is the case for the local heat equation. We are convinced that it is 
also the case for compactly supported kernels but the lower estimate of ui which implies a a strictly 
less than p does not allow us to conclude. 

(iv) The same arguments as in Remark 14.41 are valid in this case, replacing the condition on the initial 
data by ui(t) * («o)- < 00 and oj(t) * (u ) + — 00 for any t > 0. 

Now we turn to uniqueness results. In order to do so, we use the fact that, given a function / ^ 0, the 
function u(t) * f is a supersolution of (|l.lj) . We shall prove first that the constructed solution remains 
inside a suitable class of solutions, and then use the supersolution provided by u(t + 1) * / in order to 
apply the comparison principle and get uniqueness. 
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Proposition 5.6 Let J be a radially symmetric kernel, compactly supported in B p for some p > 0. Let 
uq be a locally bounded initial data, a G (0, l/p) and cq > 0. If \uq(x)\ ^ coe"' 21 ' 1 "' 21 ', there exists A > 
such that 

\u(x,t)\ ^ Ce^ ln ^. 

for any [3 £ (a, l/p). 

Proof. Let us split again the initial data into its positive and negative part, |mo| = (uo)+ + [uo)-- Then 

u+(x,t) = e-*(uo)+(a!) + (w(t) * (uo)+)(a:) < e "'co e a ^ + c&w(t) * e a l x l ln W 

and the same holds for u~ . For £ € (0, T), we have the estimate oj[x,t) ^ ce~( 1 / p ^ x \ ln \ x \ +c ( T ^ x \ , for 
some c > and c(T) e K. Hence, 

(w(t)* e a l x l ln l x l)(a;) sC c' / e - (1/p)|y|ln|!/l+c ' (T)l!y|+a|;i; ^ |ln|;i: ^ 1 dy, 

^ c ' / e -(l/p)|*|Jn|K|+c'(T)|K|+a|x|la|z|+Qi|»Hn|y|+a|x|+a|l/| (Jy 

^ ^galzllnlal+alxl /" e (a-l/p)|y| ln |i/|+(c'(T)+a)|»| dy ^ 
< q" q ot\x\Xn |x|+a|x| 

Summing up, we get that for all (3 £ (a, l/p) and x large enough, 

u + (x,t) < C , e a l a: l ln l x l+ a l x l < e^H^I and u~{x,t) < Ce a|a:| ln N+ Q N <; e ^ x IH*l 
Finally, we note that as T is arbitrary, 

\u{x,t)\ < Ce* |ln|11 

for all i > 0. 



□ 



In the following, we denote by fp the function (x) = e ' 



'\x\ ln \x\ 



Lemma 5.7 £ei J be a radially symmetric kernel, compactly supported in B p for some p > 0. Let uq 
be a locally bounded function, a £ (0,1/p) and Co > 0. If \uq(x)\ $J Co e a W ln l I l ; then for any T > 0, 
/3 £ (a, l/p) and A > Zarge enough the function 

rl>(x,t) :=A(Gj(t + l)*f )(x) 

is a supersolution of ([L~T ]) - (fL~2 ]l on R N x [0,T]. 

Proof. We assume here that p = 1 in order to simplify the proof a little bit, the modifications being 
straightforward in the other cases. Since uj satisfies (|2.3[) . the function ip verifies 

ipt = J * ip - ip + Ae ~ (t+1) J *ff3^J*?p-ip, 

so that we have a supersolution of the equation, and we need now to compare the initial data. 
We start from the fact that for t £ (0, T), and any a > 1, we have the estimate 

w(aj, t + 1) > ce- (lM|11 ln N+^W , 
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for some c > and c(T) € K. Hence, 

(w(t+l)*/ /J )(a;)>c / e -<V*)lvlMi/l+c(T)|v|+0|«-y|M*-y| dl , j 



^ c / e -(iMIal In |2/|+c(T)|y|+/3| K | In |x| dy ^ 

J{\x-y\>\x\} 

^ c(a:)e^'|i"W , 

where the function c{x) is given by 

C{x) = cf e -{l/°)\v\ln\v\+c{T)\y\ dy _ 
J y|>|x|} 

We claim that the constant c(x) is uniformly bounded from below. Indeed, since 

\x-y\ 2 = \x\ 2 + \y\ i -2^2x i y i , 

the set {\x — y\ > \x\} contains at least {yiXi < for i = 1, • • • ,N}. Moreover, since the integrand is 
radial and integrable, we have 

C(x)>c[ e -Wn\y\+c(T)\y\ dy = c f e -Wn\v\+c(T)\y\ dy = Ci >Q ^ 

J{Vi, x iyi <0} J{\/i, yi>Q} 

and thus we get that 

(u;(x,t+l)*f )(x)>c 1 e^ bl \ x \. (5.4) 
This implies that if A > cq/ci , then ip(0) — (w(l) * fp) ^ u . □ 

Theorem 5.8 Let J be a radially symmetric kernel, compactly supported in B p for some p > 0. Then 
for any < a < 1/p, there exists a unique solution of fj 1 . 1 1) (ll.2[) such that \u(x,t)\ ^ C(t) e"' 1 ' ln l x l , 
where C(-) is locally bounded in [0,oo). 

Proof. Notice first that by assumption, the initial data satifies |tto| ^ C(0) e ln l x l so that we know 
how to construct a solution satisfying the assumptions of the theorem. Now, let us fix f3 G (a, 1/p) and 
consider the supersolution 

i>(x,t):=A(u>(t + l)*fp)(x), 

which grows at least as e ^'^ ln ' x ' on x [0, T] - see (15.41) . Since this supersolution grows strictly faster 
than any solution satisfying t)| ^ C(t) e a l x ln \ x < , we can apply Theorem 13.21 to get uniqueness. □ 

5.2 Gaussian kernels 

We assume here that the J is a centered gaussian with variance a 2 = l/(27 2 ) (we write it under this 
form in order to simplify the statements), that is: 

J(y) = c( 7 )e-T |y| , c( 7 ) = ( jf e-T M dy) = (-L) . (5.5) 

We use the same method as in the case of compactly supported kernels to estimate the nonlocal heat 
kernel associated to the equation. Let us first give a slightly better estimate than that of Lemma [ 
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Lemma 5.9 There exist c > and fj, G (0, 1) depending only on J such that 

Vn>l, VxeB n , J* n {x) > c/i". 

Proof. The method is exactly the same as for Lemma liOl except that since J is not compactly supported, 
it is even easier. In the present case we set 



fi:= / J(y)dy>0, 

J{3/2<y N <2} 

and it is enough to check that if |jc| < n, and 3/2 < < 2, then \x — y\ 2 ^ (n — 1/2) 2 + 2 . Hence for n 
big enough, we have for such x, y 7 that \x — y\ ^ n and we can proceed by induction as in Lemma l5.3l □ 

Proposition 5.10 Let J be defined by (I5.5|) for some 7 > 0. Then there exist some positive constants 
Ci, C2, C3, C4 > depending only on 7 and N such that for \x\ > 1, 

cie-'e " 7|a;|(ln M)1/2 e ( C2+ln '^l < u(x, t)^c 3 e l x D 1/2 e (<*+ ln 01*1 . 

Proof. Fixing first |x| > 1, we begin by using the explicit formula for convolution of gaussian laws: 

x N/2 2 3 

J* n ( y ) = c( 7 )«(J_) e ^ l«l /" 
where 0(7) is defined in (|5.5[) . Then it follows that we can estimate from above 

K 



where JsT > 1 is an integer to be fixed below, depending on a;. Using (j5.2[) for the second term, we get 
that for some constant C(7, AT) > we have: 



w(z,t)<C( 7j J\0 (c(7)t)^(e 



-7 2 |a;| 2 /^ 4. J_ 

if! 



Now, in order to optimize this estimate as |x| — > 00, we choose K such that both terms are comparable, 
which means as |x| — > 00, using Stirling's formula: 

j 2 \x\ 2 ~ K\n(K!) ~ K 2 \nK. (5.6) 

Taking logs, we get In \x\ ~ In if and then taking the square root in (I5.6[) we obtain 

XlnA:-7|a;|(ln|a;|) 1/2 . 

With this choice, the estimate becomes: 

u(x, t) < 2C( 7 , iV) (c( 7 )i) K e * ^ < 2C( 7 , A 7 ) (c( 7 )i)* e " K ln *+° ^ , 

and we argue as in Lemma 15.21 Notice that the coefficients C2 and C4 include the constant 0(7) since 
\n(c(j)t) = lnc(7) + Int. 

To get the lower bound, we first notice that for any fixed n, 



uj(x, t) ^ e 



_ t t n J* n (x) 
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Then we choose n = [K(x)], the entire part of K defined by (|5.6I) . we use Lemma [531 

and we end up as in Proposition l5.il □ 
As a direct consequence we have: 



Theorem 5.11 Let J be defined in (|5.5p for some 7 > 0, let uq be a locally bounded initial data and 
Co > arbitrary. Then the following holds (we consider only \x\ > 1): 

(i) If \uq(x)\ ^ Co e a l x K ln 1 for some a < 7, then there exists a global solution of (jTTTJ) . 

(ii) Ifu (x) > CQe m(^\) 1/2 for 

some /3 > 7, then there does not exist any nonnegative solution of (jl.ll) . 

(Hi) If uq(x) = Co e 7 ' 31 "' 111 ' 1 !) 1 where for some a < < (3, as f(s) ^ fis, then the minimal 

solution of blows up in finite time. 

Proof. Points (i) and (ii) are done exactly as in the case of compactly supported kernels. It remains to 
show (Hi). Since as ^ f(s) ^ /3s, then we know that at least, for t > small enough, the convolution 
ui(t) * uq converges. More precisely, it is enough to choose /3 + C2 + lni < 0, hence t < e~^~ C4 . On the 
other hand, if t is big enough, using the estimate from below of oj(t), we know that the convolution blows 
up, at least for t > e _Q_C 2 , Hence the solution is defined for a short time interval, but it eventually 
blows up in finite time. □ 

Concerning changing sign solutions and uniqueness, the results for the compactly supported case have 
a direct translation in the present situation. The adaptations being straightforward, which essentially 
consist in changing ln(|x|) to ln(lxl) 1 / 2 for |x| > 1, we skip the details. 

In order to avoid the problem with the log if |x| < 1, we consider the function 

g a (x) := e«l*l(MM+i)) 1/2 . 

Everything follows from the estimates in Lemma 15.71 which are valid in the following form here: 
Lemma 5.12 Let a € (0,7) and let Uq be a nonnegative continuous function such that 

uo(x) <coe a N^ 1+ W» 1/2 . 

Then for any T > 0, if A > is big enough the function tf)(x, t) :— A(u(t + 1) * g a )(x) is a supersolution 
of (fTTTT) - (fTT§|) on R N x [0,T]. Moreover, for t € [0,T], there exists a constant A' > such that 

1>{x,t) ^' e «M(Ml+M)) 1/2 +«l*l. 

Proof. The proof follows exactly that of Lemma 15. 71 the only modification we need is provided by the 
inequality (a + b) 1 ^ 2 < a 1 / 2 + b x l 2 for the upper estimate, applied to (ln(l + |x|) + ln(l + lyl)) 1 / 2 . □ 

With this Lemma, Theorem 15.81 can be written here as follows, with obvious adaptations: 

Theorem 5.13 Let J be defined in (|5.5[) for some 7 > 0. Then for any < a < 7, there exists a unique 
solution of (ll.lD - (|1.2j) such that \u(x,t)\ ^ C(t) e a \ x \ ln l x l , where C(-) is locally bounded in [0, 00). 
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6 Explicit solutions. Asymptotic behaviour 



Let J be a kernel such that 

70 := sup {7 > : J J(y)M 7 dy < 00} e (0, 00) (6.1) 

Since J is symmetric it is easy to see that for any integer p € (0,70) the following computation makes 
sense: 

J * \x\ 2p - \x\ 2p = [ J(y)(\x - y\ 2p - \x\ 2p ) dy 



E( 2 2 -)N 2l / J(2/)|y|2(p " l) d2/ = £ 

i — n V / <* i— n 



i=0 x ' " i=0 

With this expression in mind we prove the following theorem: 



m 2 ( p - l )\x\ 21 . 



Theorem 6.1 Under hypothesis (|6.1[) . for any integer p £ (0,7o) the unique solution of (II. lj) with initial 
data uq(x) = \x\ 2p has the explicit form 



i(x,t)^\x\ 2 P + j2^(x)j, (6.2) 

k=l 



where 



(x) = J*\x\ 2p -\x\ 2 r, Ck (x) = ——(j*c k - X {x)-c k ^{x)). (6.3) 

(It i) 

Proof. We first check that u given by (|6.2[) is a solution of ([l.ip if the coefficients c k are given recursively 
by ([O]). Indeed, 



Ut ='^^Ck(x)t k ~ 1 , and J * u — u = J * uq — uq + — ( J * c k {x) — c k {x)), 



k=l fc=l 

form where we get that u satisfies the equation by replacing the coefficients by its recursive expression. 
Notice also that u(x,0) = uq(x). □ 

Examples. Consider the initial data uo = \x\ 2 . The solution u of (jl.ip has the explicit form 

u(x, t) — \x\ 2 + rri2t. 

For Mo = |x| we have 

u(x, t) = \x\ 4 + ^mjt 2 + (m 4 + 2m 3 \x\ 2 )t. 

Remark 6.2 It follows form (JO]) that u(x,t) = \x\ 2p + t p + o(t p ) as t ->• 00 which is the same relation 
between space and time as for the heat equation. 

7 Extensions and comments 
7.1 Back to the initial trace 

There is a huge difference between the nonlocal case treated here and the heat equation. Strong non- 
negative solutions have an initial trace which is necessarily a locally integrable function, not a general 
measure. This has to do with the absence of regularization of this nonlocal equation. 
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However it is not clear what happens for changing sign solutions, since the monotonicity property 
that was used in Proposition 12.51 is not valid here. In the case of the usual heat equation, it is possible 
to construct a dipole solution with initial data uq = (So)'. A similar construction can be done here, but 
the associated solution maintains the derivative of the delta measure, as can be easily seen passing to 
the limit in the explicit formula for the solution. This is a hint that probably the initial trace of strong 
changing sign solutions is also a locally integrable function, but we are not able to prove it so far. 



7.2 Singular kernels 

In [T], the authors address the case of the fractional Laplacian, J(z) — l/\z\ N+a , from the point of view of 
unbounded solutions. To be precise, they consider the more general equation u t +H(t, x, u, Du)+g[u] = 0. 
The Hamiltonian H has to satisfy some structure assumption and the nonlocal term g[u] is given by 

g[u){x) :=- f {u(x + z) - u(x) - Du ^ ' Z } d f i(z) , (7.1) 

J R N I 1 + \Z\ 2 J 

where /i is a symmetric Levy measure with density J satisfying: J inf (|z| 2 , V)J(z) dz < oo. Of course this 
equation has to be understood in the sense of viscosity solutions. 

We already pointed out that, as far as finding the optimal behaviour of the initial data is concerned, 
only the tail of J plays an important role, not its singularity at the origin. This is reflected in jj] in the 
fact that the integrability condition only concerns the bounded part /if, of the Levy measure (defined as 
the restriction of /i to the complement of the unit ball). Thus, our various estimates for the optimal data 
remain valid even for singular measures. 

For instance, if we consider the equation 

Ut(x, t) = / {u(x + z) — u(x)} J(z) dz , 



where J may be singular at the origin and behaves like a power for large \z\, then existence and uniqueness 
of solutions is based on the existence of a function / such that 



Af(z) := / \f(x + z)- f(x)\j(z) dz < A/(:r) 

.TON I. J 



This is what is done in the construction of barriers in [T] . Here is also another explicit situation: following 
Section H31 let the kernel J be a tempered a-stable law with < a < 2, often used in finance modeling: 

J ( x ) = rzjwTn ■ 



Then for any 7 < 1, the function f(x) = e 7 ' 2 ' is a barrier. Indeed, 



it/ C e 7 l z l - 1 f e 7|z| - 1 

^^L*>W e ds + i>«>TF^~'" lds 

< e^(C 1 +C 2 ) = Xf(x). 

Notice that to give sense to those integrals near z = 0, either one has to take a < 1, or to understand 
the operator as the principal value of the integral if 1 ^ a < 2. 



l\x+z\ _ e y\x\ r e 7|x+z|-7|x| _ ]_ 
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A similar computation shows that for any < a' < a, the function /(x) = e'^l + \x\) a is also a 
barrier. Then, existence and uniqueness results that complement those of [I] can be obtained provided 
the initial data satisfies |uo(aO| ^ e 1*1(1 + |x|) Q for some < a' < a. 

In the introduction of pQ, it is conjectured that it may always be possible to construct some barrier for 
any given Levy measure. This may be true, however we point out again that in the case of fast decaying 
(or compactly supported) kernels, the elliptic barriers do not give the optimal behaviour. 

Typically, if we think about the singular measure with density 

!{M<i} 



H(z) 



N+a 



then the barriers that can be constructed allow to get existence only up to an exp (|x| ) -behaviour, but 
they miss the exp(|x| In |x|)-behaviour. 

We finish this section by mentioning that a way to get optimal existence for singular measures is to 
to approach the singular measure with a monotone sequence of bounded measures. Passing to the limit 
has to be done in the viscosity sense of course, and this allows to get existence of a solution for a general 
Lvy measure \i provided the initial data satisfies the estimates for the nonsingular part of \i. 



7.3 Large deviations 

In [5] and [3] , some bounds for the same nonlocal linear equation were derived, but related to a different 
problem: estimating the error when approaching the solution u in all M. N from solutions ur of the Dirichlet 
problem in a ball -Br, as R — > oo. 

The two problems are somehow related in the sense that measuring the difference (u — ur)(T) amounts 
to measuring the total amount of sample paths that can escape the ball Br between times and T. And 
in some sense, this is another way to estimate the nonlocal heat kernel associated to the equation. 

Thus, we recover for instance the typical (e ~- Rln -^-bounds for compactly supported kernels and 
( e _i? )-bounds for exponentially decaying kernels. In the case when J decays at most exponentially, more 
bounds are to be found in [3] which should give the behaviour of the initial trace for various kernels 
J, according to their decay. 
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